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Abstract 

The electromagnetic potential consisting in the Coulomb plus the magnetic 
moment interactions between two nucleons is studied in nucleon-deuteron 
scattering. For states in which the relative N — d angular momentum L has 
low values the three-nucleon problem has been solved using the correlated 
hyperspherical harmonic expansion basis. For states in which the angular 
momentum L has large values, explicit formulae for the nucleon-deuteron 
magnetic moment interaction are derived and used to calculate the corre- 
sponding T-matrices in Born approximation. Then, the transition matrices 
describing N — d elastic scattering have been derived including an infinite 
number of partial waves as required by the behavior of the magnetic 
moment interaction. Appreciable effects are observed in the vector analyz- 
ing powers at low energies. The evolution of these effects by increasing the 
collision energy is examined. 

I. INTRODUCTION 

The study of the magnetic moment interaction (MM) in the two-nucleon (2A^) system 
has been subject of many investigations (see Refs. [1,2] and references there in). Although 
the intensity of this interaction is very small compared to the nuclear interaction, its long 
range behavior produces significant effects in nucleon-nucleon (NN) scattering. Almost 
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all modern NN potentials have been constructed considering the electromagnetic (EM) 
interaction used in the Nijmegen partial-wave analysis which includes the MM interaction 
between the two spin-i particles as well as corrections to the p — p Coulomb potential as 
two-photon exchange, Darwin-Foldy and vacuum polarization terms. When 2N scattering 
observables are computed with one of these potentials the long range behavior of the EM 
interaction implies an infinite sum in the partial-wave series. For the particular case of the 
MM interaction, in Refs. [1,2] it has been shown how to sum analytically these infinite series 
for p — p and n — p scattering. Important effects of the MM interaction has been observed 
in both n — p and p — p vector analyzing powers at low energies. 

Due to the fact that 2A'" potentials are constructed by fitting the NN available data, 
the three-nucleon {3N) system is the simplest one in which these potentials can be used to 
make predictions. However, in the description of the 3N continuum the MM interaction and 
corrections to the Coulomb potential has been systematically disregarded. This omission 
has been justified in the past by the intrinsic difficulties in solving the nuclear problem. At 
present, the 3N continuum is routinely solved by different techniques making possible the 
treatment of those electromagnetic terms beyond the Coulomb interaction. 

In the present paper we study N — d elastic scattering including Coulomb plus MM in- 
teractions. Previous description of this process without considering the MM interaction has 
been performed by the authors using a technique based on the Kohn variational principle 
(KVP) [3,4] and expanding the scattering wave function in terms of the correlated hyper- 
spherical harmonics basis [5,6]. Following these works we perform a partial- wave decompo- 
sition of the scattering process. For states with low values of the relative orbital angular 
momentum L of the projectile and the target, the process is studied by solving the complete 
3A^ problem with the Hamiltonian of the system containing nuclear plus Coulomb plus MM 
interactions. For states with L values sufficiently high, the centrifugal barrier prevents a 
close approach of the projectile to the target. So, the collision can be considered peripheral 
and treated as a two-body process. Furthermore, in these states only the EM interaction 
gives appreciable effects and the corresponding scattering amplitudes can be calculated in 
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Born approximation. The value of L at which the treatment of the problem changes from a 
three-body description to a two-body description is to some extent arbitrary and could be 
different at different energies. In practice it can be taken equal to the maximum L value 
considered when the problem is solved neglecting the MM interaction. 

We apply this procedure to calculate the ?>N vector analyzing powers where the main 
effects of the MM interaction can be observed. For p — d scattering a sizable increase in Ay 
and iTii has been obtained at low energies which is, however, insufficient to explain the usual 
underestimation produced by modern A^A^ forces [7,8]. Other observables as the differential 
cross section and the tensor analyzing powers suffer minor modifications, of the order of 1% 
or less. For n — d scattering a pronounced effect at very small angles is observed. In fact, 
the scattering amplitude has a term sin^/(l — cos^) which diverges for ^ — > similarly to 
the n — p case [2]. The experimental observation of this divergence is problematic since it 
occurs at extreme forward angles (a fraction of degree). This is different from the p — d case 
in which the Coulomb divergence dominates in that region. Regarding the vector analyzing 
powers, the MM interaction tends to slightly flatten the n — d Ay around the peak and to 
produce a pronounced dip structure at small scattering angles. 

The importance of the EM interaction in the description oi N — d scattering decreases 
as the energy of the process increases. Around Eiab = 16 MeV the improvement given by 
the MM interaction at the peak of Ay and iTu ior p — d scattering is already less than 5%. 
On the other hand Coulomb effects are important below Eiab — 30 MeV [3] . Here we show 
that at Eiab = 65 MeV they are considerably reduced in most of the observables with the 
exception of T21 where still some effects can be observed. This analysis will serve to justify 
the application of standard n — d calculations to the description oi p — d scattering at high 
energies [9]. 

The paper is organized as follows. In Section II the N — d MM interaction is given. The 
corresponding T-matrices are calculated in Born approximation for both n — d and p — d 
scattering and final forms for the transition matrices are given. In Section III the transition 
from a 2)N description to a description is discussed. It is shown that the 3A^ T-matrix 
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tends to the 2A'" T-matrix as the value of L increases. In Section IV the vector analyzing 
powers are calculated including the MM interaction and compared to the available data. 
The differences between the theory and the experiments around the peak of the observables 
are analyzed. In Section V we present our conclusions. In the Appendix the n — d as well 
as the p — d MM interactions as two distinctive particles are derived. 

II. MAGNETIC MOMENT INTERACTION 

Following the notation used in the determination of the Argonne vi^ (AVI 8) poten- 
tial [10], all modern NN potentials can be put in the general form 



The short range part v^{NN) of these interactions includes a certain number of pa- 
rameters (around 40), which are determined by a fitting procedure to the A'^A'^ scattering 
data and the deuteron binding energy (BE), whereas the long range part reduces to the 
one-pion-exchange potential v'^{NN) and the electromagnetic potential v^^{NN). 

The AV18 potential includes the same v^^{NN) used in the Nijmegen partial- wave 
analysis except for short-range terms and finite size corrections. The v^^{pp) consists of 
the one- and two-photon Coulomb terms plus the Darwin-Foldy term, vacuum polarization 
and MM interactions. The v^^{np) interaction includes a Coulomb term due to the neutron 
charge distribution in addition to the MM interaction. Finally, v^^{nn) is given by the 
MM interaction only. All these terms take into account the finite size of the nucleon charge 
distributions. Explicitly the two-nucleon magnetic moment interaction in the center of mass 
reference frame reads: 



v{NN) 



EM 



(NN) + v^'iNN) + v^{NN) . 



(1) 



VMM{np) 



vmm{pp) 




(2) 
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2MnM, 



a 




(L • S + L • A) 



(3) 



VMM{nn) 




a 




(4) 



In the above formula F5, Ft and describe the finite size of the nucleon charge distri- 
butions. As r — > 00, ^ whereas Ft ^ 1 and F/^ — > 1. Mp (M„) is the proton (neutron) 
mass and M„p is the n — p reduced mass. The MM interaction presents the usual be- 
havior and has an operatorial structure with a spin-spin, a tensor and a spin-orbit term. In 
the n — p case, this last term includes an asymmetric force (proportional to A ={cri—(Tj)/2) 
which mixes spin-singlet and spin-triplet states. This term is expected to have a very small 
effect. 

The EM interaction has been studied in the description of bound states in ^4 < 8 nucleon 
systems [11]. Recently a detailed analysis of the contribution of the electromagnetic terms 
to the ^He-^H mass difference has been performed [12]. A first analysis in three-nucleon 
scattering has been done by Stoks [13] including the MM interaction in Born approximation 
at high L values. However, the T-matrices used at low L values were calculated without 
considering the MM interaction. In this approximate treatment of the process the main 
modifications were obtained in the n — d vector analyzing powers at forward angles. No 
modifications were observed in other observables as the differential cross section and tensor 
analyzing powers and in the maximum of Ay and iTii. As a consequence, the conclusion 
was that the MM interaction does not improve the theoretical underestimation of the last 
two observables. However, disregarding the MM interaction could not be correct in the 
description of low partial waves which govern the polarization observables at low energies. 
In Refs. [14,15] the MM interaction has been included in the calculation ol N — d scattering 
observables. However in these analyses its contribution was limited to a low number of 
partial waves. The contribution from waves with high L values was neglected. In the 
present paper we will include the MM interaction in both regimes in order to perform a 
complete description of the collision process. 

For the case A — 2, the contribution of the MM interaction to the scattering amplitude 
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has been extensively studied [1,2]. It has been shown that due to its behavior the 
scattering amphtude results in a slow convergent series whose leading term can be summed 
analytically. A similar analysis can be performed for N — d scattering. The starting point is 
the A'" — d transition matrix M which can be decomposed as a sum of the Coulomb amplitude 
fc plus a nuclear term, namely 

Mi^:{0)^f,{e)5ss'S.u' + ^ E V2LTliL0Su\Jiy){L'M'S'u'\Ju) 

L,L',J 

exp[i{aL + av - 2a^)] 'tH', Yvm'{Q. 0) . (5) 

This is a 6 X 6 matrix corresponding to the two possible couplings of the spin 1 of the 
deuteron and the spin 1/2 of the third particle to S^S' — 1/2 or 3/2 and their projections 
V and v' . The quantum numbers L, L! represent the relative orbital angular momentum 
between the deuteron and the third particle and J is the total angular momentum of the 
three-nucleon scattering state. "^Tff,' are the T-matrix elements corresponding to a Hamil- 
tonian containing nuclear plus Coulomb plus MM interactions and ui, are the Coulomb 
phase-shifts. The n — d case is recovered putting fc = cfl = 0. When the MM interac- 
tion is not considered the sums over L, L', J converge very fast due to the finite range of 
the nuclear interactions. Typically in the low energy region {Eiai, < 50 MeV) states with 
L, L' > 10 can be safely neglected. However, when the MM interaction is considered, an 
infinite number of terms contributes to the construction of the scattering amplitude. In this 
case the sums on L, L' can be divided in two parts. For L, L' < L^ax the T-matrix elements 
correspond to, and are obtained from, a complete three-body description of the system. For 
L, L' > Ljnax the centrifugal barrier is sufficiently high to maintain the third particle far 
from the deuteron and the description of the state can be performed as a two-body system. 
In general L^ax can be fixed in such a way that when the collision proceeds in states with 
L, V > Ljjiax the nuclear interaction can be safely neglected and only the Coulomb plus MM 
potentials contributes to the N — d scattering. It is therefore convenient to introduce the 
MM interaction between a micleon and the deuteron as distinct particles. Its specific form 
can be obtained summing the MM interaction between each nucleon of the deuteron and 
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the third nucleon at large separation distances. Alternatively, the N — d MM interaction 
can be obtained directly in one- photon exchange approximation between a spin-1 and a 
spin-1/2 particle from a non-relativistic reduction of the corresponding Feynman diagram. 
Here below the MM n — d and p — d interactions are explicitly given. The details of the 
derivation are reserved to the Appendix. 

VMM{nd) = -y,[-^SL + ^]t£l^(L • Sn. + L . A„,)] , (6) 
VMMipd) = -^[^^i + - ^)(L . S,. + L . A,,) 

Sf,, = 3{SN-r){S,-r)-SM-Sa, N^n,p (8) 
S'/ = 3(S, • f)2 - 2 , (9) 

where is the deuteron mass, Mjvd is the corresponding nucleon-deuteron reduced mass 
and /Id, Qd are the magnetic and the quadrupole moments of the deuteron, respectively. 
Moreover, Sat^ = Sat + whereas Aat^ = Sat — S(i. The deuteron- nucleon distance is r and 
r is the unitary vector giving their relative position. 



A. n — d case 



Let us first discuss n — d scattering including the MM interaction. For relative states 
verifying L, L' > L^ax the description proceeds as a two-body process and the T-matrix 
elements corresponding to a. n — d state with total angular momentum J, relative angular 
momentum L and total spin S are given in Born approximation as 

'T^i'' = -k{'^)<^L'S'j\vMMind)\nLsj> . (10) 

The relative motion of the n—d system is described by the regular free solution of Schrodinger 
equation 

^LSJ^jL{kr)[YL{f)(^Xs]jj. , (11) 
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with k"^ — {2Mnd/'h^)Ecrm 3l a spherical Bessel function and xs the total spin function. 

The T-matrix elements corresponding to the spin-orbit term of the MM interaction 
proportional to L • S + L • A are 



'''^SS' — ~^ ^so < ^L'S'j] ^ l^LSJ > 



-k Cfi 



with 



2.932 X 10-'^ fm 



. r -^ML, , (12) 

2L(L + 1) ' ^ ^ 



(13) 



and 



= (_l)i+^+s-s'-i/2^6(2^+ 1)(25' + 1) 



I S' 1 
S i 1 



^L(L + 1)(2L+1 




The T-matrix elements of Eq.(12) can be used in Eq.(5) for values of L, L' > Lmax- 
Moreover, for fixed values of L the sum over J can be performed analytically using sum- 
mation properties of Clebsh-Gordan coefficients. The convergence of the sum on L is slow 
enough to prevent a safe truncation of the series. Therefore, after summing all terms for 
L > Lmax-i the contribution of the spin-orbit term to the transition matrix of Eq.(5) results: 



sin 6 



L=i L{L + 1) 



(15) 



1 — cos 9 

Pl{cos9) is a generahzed Legendre polynomial and the following property has been used to 
derive the above equation 

- (2L+21 sin^ 

Moreover 



1 — cos 9 



(16) 



- S' 1 

K^s; ^ [-i)S-s'+-H^^{2S+l){2S' + 1) { ' 

S \ 1 



S S' 



-M 



(17) 



As a consequence of the r ^ behavior of the MM spin-orbit interaction a term propor- 
tional to sin ^/ (1 — cos 9) appears in the transition matrix. This term produces a divergence 



in the differential cross section at extreme forward angles and a pronounced dip structure 
in the vector analyzing powers. 

A similar analysis can be done for the term proportional to the tensor operator in the 
n — d MM interaction. The corresponding T-matrix elements are 



nd 



(18) 



with 



The angular-spin and radial matrices are 



1.675 X 10"'^ fm 



\ 1 S' 
i 1 s 



112 



X J30(2L + 1)(2L' + 1){2S + 1){2S' + 1) { 



V S' J 




u 


2 L'^ 


< 


> 






S L 2 










and 



III' = < 



2L(L+1) 



6(L+l)(L+2) 

5r_o Tl 



(19) 



(20) 



(21) 



l_ 6(L'+l)(L'+2) 

Again for fixed values of L and L' the sum over J in Eq. (5) can be performed analytically 
and the contribution to the transition matrix is 



M^^:{t) = -v^Ct^30(25 + 1){2S' + 1) < 



X {'2L + 1)V2L' + IIll' 




yL'M(^,0) . (22) 



Three different sums can be constructed corresponding to \M\ — 0,1,2 that can be 
summed numerically term by term. The convergence of the series is rather fast and a few 
tens of terms are sufficient. 

In conclusion, the n — d transition matrix including the nuclear plus the MM interaction 
can be constructed as a sum of three terms 

M'J{0) = V 5: E^2LTT(L05z/ I Jz/)(L'MW I J^) ^Tff,' YyM'{e,Q) 

L,L' J 

+M^S'iso)+M^S'it) . (23) 

When the MM interaction is neglected only the first term contributes to the transition 
matrix. When the MM interaction is included, the T-matrix elements for L, L' < L^ax 
are different from the previous case. In addition the last two terms in Eq.(23) have to be 
included. We stress the fact that the value of L^ax can be taken in such a way that for 
L, L' > Ljnax the nuclear interaction gives a negligible contribution to the scattering process 
and the interaction between the incident particle and the target is only electromagnetic. 
Typical values for L^ax are discussed in Sec.IV. 



B. p — d case 

As for the n — d case, the T-matrix elements corresponding to a two-body description of 
the p — d system with total angular momentum J, relative angular momentum L and total 
spin are given in Born approximation 

'Tli: = -k{^^) < nL>s'j\vMM{pd)\nLSJ > ■ (24) 
n 

Here the relative motion of the p — d system is described by 

^LSJ = FL{r], kT)[YL{f) ® xs\jj. (25) 

with k"^ — (2Mpd/h'^)Ecm, Fl a regular Coulomb function and rj the usual Coulomb param- 
eter. 



10 



Let first consider the spin-orbit terms of the MM interaction in Eq.(7) proportional to 
(L • S + L • A) and (L • S — L • A). The following matrix elements entering in the calculation 
of the T-matrix are defined 



(26) 



with [16] 
III = 



+ 



2L(L + 1) 2L(L + 1)(2L + 1) 



rjTT + 1 + rjTTCothrjTT — 2rj^ ^ 



(27) 



In Eq.(26) the angular-spin matrix '^Mgg,{+) is equal to the matrix '^Mgg, defined in 
Eq.(14), whereas 



^M|5,(-) = (-l)'^+^-V2^24(25+l)(25'+l) < 



/ 1 



1 S 
S 1 1 



y'L(L+l)(2L+l) 




Following Ref. [1] we isolate the first term of III and proceed toward a summation of 
the related amplitude as we have done for the n — d case. The corresponding contribution 
to the transition matrix of Eq.(5) for L > L^ax results 



M5 (so) = -[CtoK'J (+) + C:^Ki', (-)] 



X 



cos 9 + 2e-'^'<"-^^ - 1 (2L + 1) ^.a^.,^) m 



sin^ 



(29) 



To get this final form we have used the following analytical summation of the series [17] 



V -^fi±i(e2--Pi(cos^) = ^[cos^ + 2e-*''i"(^) _ i] , 
L[L + 1) sinb* 

which can be obtained from the series of the Coulomb amplitude 



(30) 



fc{0) = E(2^ + ^K^'""' - I)^l(cos^) = -2t7]- -e 

r n -1- — COS U 



— ir;ln(- 



(31) 



L=0 



using the recurrence relations of the Legendre polynomials and the following relation of the 
Coulomb phase-shifts 
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g2i(ai-ai_i)(^ - i77) = L + i77 . (32) 

In Eq.(29) K^^' {+) = K^^,' defined in Eq.(17) and 



K5\-) = (-irVl2(25 + l)(25' + l) 



1 5' i 
5 1 1 



^ S S' 1 ^ 
1/ -u' -M ^ 



^\M\,i ■ (33) 



Moreover 



'-'so — 


-aMpd 1 








-aMpd 


( l^d 





^ -3.775 X 10-3 fm , (34) 

^-^1 -5.936 X 10"^ fm . (35) 

Tfie term proportional to Cj^ is mucfi smaller due to the small magnetic moment of the 
deuteron. The same happens to the term proportional to Qa in Eq.(7) due to the small 
quadrupole moment of the deuteron and will not be discussed here. The analysis of the 
term proportional to the tensor operator in the MM interaction proceeds similarly to that 
one performed in the n — d case, taking care that now the radial integral III is given by 
Eq.(27) and Il,l+2 — l\L + 1 + ir]\~^\L + 2 + irj]"^ [16]. In conclusion the transition matrix 
can be constructed as a sum of different contributions 



LU J 

X exp[i((7L + av - 2(7o)] 'tI^, Yl'm'{0, 0) + M^^' M + B^J , (36) 

where M^^,' (so) is defined in Eq.(29) and B^^,' includes the contribution of the remaining 
terms in Eq.(27) and those coming from the tensor operator. The B^^' matrix elements 
can be evaluated summing the corresponding series numerically for L, L' > Lmax until 
convergence is reached. 

III. THE 3A^ AND N-D T-MATRICES IN BORN APPROXIMATION 



The calculations of the obscrvables in N — d scattering can be obtained from the tran- 
sition matrices of Eqs.(23) and (36). Accordingly, after a partial wave decomposition, it is 
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necessary to calculate the three-nucleon T-matrices for states with total angular momen- 
tum J in which the deuteron and the incident nucleon are in relative motion in the regime 
L < Lmax- As discussed before, states having L > Lmax are described as a two-body pro- 
cess. Therefore it is appropriate to make a link between the two regimes and show in which 
manner the three-nucleon T-matrix smoothly tends to a two-body T-matrix as L increases. 

The KVP in its complex form establishes that the T-matrix elements are functionals of 
the three-nucleon scattering state 



The stationarity of this functional with respect to the trial parameters in the three-nucleon 
scattering state ^lsj required to obtain the T-matrix first order solution. The second 
order estimate is obtained after replacing the first order solution in Eq.(37). In this formal- 
ism [18] the continuum state is usually written as a sum of three Faddeev-like amplitudes, 
each of which consists of two terms: 



here Xj, are the Jacobi coordinates corresponding to the i-th permutation of the particles 
indices 1,2,3. The first term, ^'c, describes the system when the three-nucleons are close 
to each other. For large interparticle separations and energies below the deuteron breakup 
threshold it goes to zero, whereas for higher energies it must reproduce a three outgoing 
particle state. The second term, ^tsj' describes the asymptotic configuration of a deuteron 
far from the third nucleon and explicitly it is: 




(37) 



(38) 



1=1,3 




'LSJ 



(x„y,) + E'^ff''^i 



(39) 



L'S' 



where 




(40) 




(41) 
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Besides a factor \/k, ^lsj function Qlsj given in Eqs.(ll) and (25) for the n — d 

and p — d system respectively, in which xs represents the deuteron wave function of spin 1 
coupled with the spin | of the third nucleon to total spin S. In fl^sj the regular relative 
function ji or Fl is replaced by the corresponding irregular solution rji or Gl regularized 
at the origin [5] . The normalization of the asymptotic states verifies 

^^j:[i^LSji^i^yi)\H-E\Qi,s,j{^j,yj)) - {nisji^i,y,)\H - E\n^,s'ji^j,yj))] = 1 , 

(42) 

M being the nucleon mass. To be noticed that in the three-nucleon process the energy in 
the center of mass reference frame is 

with the deuteron ground state energy. Moreover, the factor l/(2-\/3) in Eq.(42) is related 
to the definitions of the Jacobi coordinates in terms of the particle coordinates: 

Xj = rj - Tk 

The Born approximation of the T-matrix is obtained from Eq.(37) replacing the wave 
function ^ by the regular function Q° and putting the first order T-matrix equal to zero: 

M 



\JrrSS'\ ^ 

For a given energy a certain value Lb exists such that for L, L' > Lb the differences 
between the T-matrix elements obtained from a complete solution of the three-nucleon 
problem or from its Born approximation are extremely small. Increasing further the values 
of L and L' we arrive to the regime L,L' > L^^x in which the contribution of the NN nuclear 
potential can be neglected. Let us consider X3,y3 the Jacobi coordinates corresponding to 
the asymptotic configuration in which nucleons (1, 2) form the deuteron and nucleon 3 is 
the incident particle. The relative coordinate between the third nucleon and the center of 
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mass of the deuteron is r^^ — (•\/3/2)y3. Starting from the above Born approximation for 
the T-matrix, the following relations are verified for L, L' > Lmax- 

['T!i:]B = -^:^Y.{^lsA^^,y^)\H - E\nl,s'ji^,,y,^ (46) 
= -3^-^ E(^Lj(x. y.) I 3) + V{2, 3)\nl,s,j{^,, ys)) (47) 
^ -3^^(OLj(x3,y3)|l^(l,3) + y(2,3)|Q°,5,,(x3,y3)) (48) 

" ^ :{nlsj{^s:y3)\vMM{Nd)\nl,s'ji^3:ys)) (49) 



Ma 

The equivalence between the second and third row is in general verified for L,L' > Lb- On 



= -2k^{QLSj\vMM{Nd)\QL'S'j)) ■ (50) 



the other hand, the equivalence between the third and fourth row is verified for L, L' > Lmax- 
In fact, Lmax can be fixed as the L value at which these two rows start to be approximately 
equal. Finally, in the last step the asymptotic three-nucleon function VP has been replaced by 
the two-body function O of Eq.(ll) once the integration over X3 and the change of variables 
ys T^Nd has been performed. In conclusion, the above approximate equalities show the 
relation between the three-nucleon T-matrix of Eq.(45) and the two-body T-matrices of 
Eqs.(lO) and (24) for high L values. 

IV. N-D OBSERVABLES INCLUDING COULOMB PLUS MM INTERACTIONS 

Elastic observables for N — d scattering can be calculated using the transition matrices 
of Eqs.(23) and (36) using trace operations [19]. The calculations presented here have been 
performed using the KVP after an expansion of the three-nucleon scattering wave function in 
terms of the pair correlated hyperspherical harmonic (PHH) basis [5,6]. As A^A^ interaction 
we have used the nuclear part of the AV18 potential plus the Coulomb and MM interactions 
defined in Eqs.(2)-(4). The asymmetric force L • A in the VMuinp) interaction has been 
included as well as the p — p Darwin- Foldy and the n — p short-range Coulomb terms. 

At energies below the deuteron breakup threshold the contribution of the MM interaction 
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is expected to be appreciable. Recently the n — d analyzing power Ay has been measured at 
Eiab — 1.2 and 1.9 MeV [20]. At these very low energies the nuclear part of the transition 
matrix (first term of Eq.(23)) converges already for L^ax = 3. The corresponding theoretical 
curves obtained using the AV18 potential, and neglecting the MM interaction, are showed 
in Fig.l (solid line). As it can be seen, the observable is not reproduced by a large amount 
which is a common feature of all modern NN forces. When the MM interaction is taken 
into account up to L^ax — 3, the analyzing powers are given by the dashed curves. There is 
a very small influence of the MM interaction in the peak of Ay with the tendency of slightly 
flattening the observable. However, this is an incomplete calculation since the inclusion 
of the MM interaction requires an inflnite number of partial waves in the calculation of 
the transition matrix. When all three terms of Eq. (23) are considered the observables are 
given by the dashed-dotted curves. It is interesting to notice the forward-angle dip structure 
which already appears in n — p scattering [2] . Only after summing the series up to oo this 
particular behavior can be reproduced. We can conclude that the MM interaction produces 
a pronounced modiflcation of Ay at forward angles but has a very small effect around the 
peak. 

In order to show the importance of the MM moment interaction in the calculations of Ay 
as the energy increases, in Fig. 2 the results at Eiab = 6.5 MeV are given. At this particular 
energy Ay has been measured in an extended angular range including forward angles [21]. 
The sohd line corresponds to a standard AV18 calculation neglecting the MM interaction and 
including partial waves up to L^ax = 8- The dashed-dotted line corresponds to a calculation 
using the AV18-I-MM potential and considering the complete series. We can observe that 
the effect of the MM interaction on the peak is practically neghgible. Conversely, it is of 
great importance at forward angles in order to describe the zero crossing. 

Besides the neutron analyzing power and the deuteron analyzing power which present 
similar characteristics, other elastic n — d observables as the tensor analyzing powers suffer 
only minor modiflcations when the MM interaction is included. The differences are of the 
order of 1% or less and they are not presented here. However when comparisons with precise 
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experimental data are performed these differences could be relevant and the MM interaction 
should be taking into account. 

For p — d scattering high precision data exist at low energies [22-25] for differential cross 
section and vector and tensor analyzing powers. Detailed comparisons to these data has 
been performed in Refs. [7,22,23,26] using AV18 with and without the inclusion of three- 
nucleon forces. In those studies the Coulomb interaction was included whereas the MM 
interaction was not. In order to evaluate the effects of the MM interaction on the vector 
analyzing powers in presence of the Coulomb field, in Fig. 3 the results of the calculations 
at Ep = 1 and 3 MeV are shown. Three different calculations have been performed at both 
energies. The solid line corresponds to the AV18 prediction neglecting the MM interaction. 
Accordingly, the transition matrix has been calculated with the first two terms of Eq.(36). 
The partial-wave series of the second term has been summed up to L^ax — 4 {Ep = 1 
MeV) and L^ax — 6 {Ep — 3 MeV). The dashed line corresponds to the same calculation 
as before but the T-matrix elements has been calculated using the AVIS+MM potential. 
The dashed-dotted line corresponds to the complete calculation including also the last two 
terms of Eq.(36). We see that the major effect of the MM interaction is obtained around the 
peak and is appreciable at both energies. There is also an improvement in the description 
of the observable at forward angles, in particular for iTu a.t Ep = 3 MeV. The observed 
modifications are due to the interference between the Coulomb and the nuclear plus the 
MM interaction and not to higher order terms, as in the n — d case, since, except for Ay 
a.t Ep = 1 MeV, the dashed and dashed-dotted line practically overlap. In fact, high order 
terms are dominated by the Coulomb interaction and the MM interaction gives a very small 
contribution. 

As the energy increases, the effect of the MM interaction on Ay and iTu diminishes as 
it is shown in Fig.4 at Ep — b and 10 MeV. Here the AV18 prediction (solid fine) has to be 
compared to the AV18-I-MM prediction (dashed line) calculated using the first two terms of 
Eq.(36) with Lmax = 8- When the last two terms of Eq.(36) are also included, the results 
are extremely close to the previous ones. As for the n — d, the tensor analyzing powers 
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present very small modifications when the MM interaction is taking into account and are 
not presented here. 

The MM interaction has difi^erent efi^ects in n — d or p — d vector analyzing powers. One 
reason is the different sign between the neutron and proton magnetic moment. Another 
reason is the interference with the Coulomb field. However the MM interaction does not 
help for a better description of the neutron Ay. On the contrary there is an appreciable 
improvement in the proton Ay as well as in iTu, in particular at very low energies. Hence 
we can examine the differences between the experimental data and the theory at the peak 
in order to see if the inclusion of the MM interaction helps to clarify a different behavior 
observed ior n — d and p — d. In Fig. 5 the relative difference [Ay{exp) — Ay{th)]/ Ay{exp) 
at the peak for n — d and ior p — d scattering is shown. In this last case both, the AV18 
and AV18+MM results have been reported. For n — d both results are extremely close at 
the peak, so the difference does not depend on which calculation (AV18 or AV18+MM) 
is considered. Without the inclusion of the MM interaction the underestimation of the 
proton Ay is much more pronounced than the neutron Ay. When the MM interaction is 
considered the difference between theory and experiment for both, p — d and n — d scattering 
are of similar size, around 25%, for all the energy values below 16 MeV. Above 16 MeV 
the differences at the peak between theory and experiment diminish. As shown in Fig. 5, 
at 18 MeV the difference is around 20%. In Fig.6 the deuteron analyzing power iTu is 
examined. The relative difference [iTii(exp) — iTu(th)]/iTu(exp) is shown at the peak for 
p — d scattering (there is no data for the n — d case) using AV18 and AV18+MM. Besides 
the first point at Eiab — 0.650 MeV which corresponds to a very small value of iTu [22], the 
underestimation of the observable oscillates around 24%, very close to the Ay case. 

Finally we wish to discuss the importance of the Coulomb effects as the energy increases. 
In fact, up to Eiab — 30 MeV we can observe appreciable differences in the description 
oi n — d and p — d elastic scattering that however tend to diminish [3]. Experimental 
data are not always conclusive since experiments with neutrons have larger uncertainties 
than those performed with protons. On the other hand, n — d calculations have been 
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often used to describe p — d scattering, in particular at high energies [9]. In order to 
clarify this approximation, in Figs. 7-8 n — d and p — d calculations at Ei^b = 65 MeV are 
compared. To make contact with the results given in Ref. [9] in which n — d scattering has 
been analyzed at this particular energy, we have consider also the Urbana IX (UR) three- 
nucleon interaction [31]. In Fig. 7 the differential cross section and Ay are shown. Three 
curves are displayed corresponding to p — d AV18 (solid line), n — d AVIS (dashed line) and 
p — d AV18+UR (dotted hne) and compared to the experimental data. In Fig. 8 the same 
calculations are shown for iTn and the three tensor analyzing powers T2o,T2i and T22. As 
expected. Coulomb effects are small at this energy. We can observe appreciable Coulomb 
effects only in T21 whereas three-nucleon interaction effects are found in the minimum of the 
differential cross section and in T21 and T22 as well. These results justify to some extent the 
description of p — d data using n — d calculations at intermediate energies, however with 
some caution in the description of particular observables. 



V. CONCLUSIONS 

The MM interaction has been included in the description oi N — d scattering at low ener- 
gies. Though its strength is small compared to the nuclear interaction, it has a very long tail 
which behaves as As a consequence, the construction of the scattering amplitude ne- 

cessitates an infinite number of partial waves. Analytical summations of the corresponding 
p — d and n — d series have been given following previous works for NN scattering. Accord- 
ingly, the 6x6 transition matrix M has been written as a sum of the standard Coulomb 
amplitude plus the MM amplitude and a finite series of T-matrices. These matrices have 
been calculated from a complete three-body description of the process with a Hamiltonian 
including the nuclear plus Coulomb plus MM interaction. For high L values, the MM am- 
plitude has been calculated as a two-body process. To this aim the MM interaction between 
a nucleon and the deuteron as distinct particles has been derived. 

Different 3A^ elastic observables have been calculated and compared to previous calcula- 
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tions in which the MM interaction was neglected. The main effect has been observed in the 
vector analyzing powers. However the modifications produced by the MM interaction do not 
improve the description of the neutron Ay around the peak. Conversely, there is an appre- 
ciable improvement in the proton Ay and in iTn, in particular at low energies. Due to the 
different effect that the MM interaction produces inn — d and p — d scattering, the relative 
difference between the calculated and the measured Ay at the peak results largely charge 
independent and approximately constant below 16 MeV. The underestimation is about 25%. 
Above this energy the difference starts to diminish. At Eiab — 18 MeV it has been reduced 
to 20% and above 30 MeV there is a much better description of Ay and iTn. This is shown 
by the calculations performed at Eiab — 65 MeV. Furthermore, we have shown that at this 
energy Coulomb effects are not important. Only T21 still shows some sensitivity. 

The main aim of this work is to describe the three-nucleon continuum using the same 
v^^{NN) used in the description of the 2A^ scattering states. In the past the MM interaction 
has been systematically neglected in the calculation of 3A^ scattering observables with few 
exceptions. Here we show how to include it and which terms are important. From the 
present analysis it can be concluded that the approximate treatment of Ref. [13] is justified 
ior n — d scattering but not for the p — d case. In fact, in the calculation of the n — d Ay the 
symmetric spin-orbit term in VMuinp) tends to depress the observable at the peak whereas 
the asymmetric term almost cancel this effect. Therefore the transition matrix of Eq.(23) 
can be constructed with the MM amplitudes M^^,' (so) and M^^,' (t) but neglecting the MM 
interaction in the calculation of the T-matrix elements '^T££,' for L, L' < L^^ax- In addition, 
the amplitude M^^,' (t) gives an extremely low contribution and can be neglected too. In 
the p — d case the interference between the Coulomb, MM and nuclear interactions does not 
allow for the omission of the MM interaction in the calculation of the T-matrix elements. 
Otherwise the improvement at low energies on the peak of Ay and iTu is lost. However, 
in the construction of the transition matrix the last two terms in Eq.(36) give very small 
contributions and, except at extremely low energies, can be omitted. 

Other small terms in the v^^{pp) interaction as the two- photon Coulomb and vacuum 
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polarization interactions have been neglected in the present analysis. These terms have 
improved the description of p — p scattering at low energies and, therefore, their inclusion 
in the description of p — d scattering is of interest. The analysis of these terms as well as 
the study of the MM interaction in p— ^He scattering is at present underway. 
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APPENDIX A 

In this appendix we briefly outline the derivation of the N — d MM interactions of Eqs. (6) 
and (9). 

We consider two particles, the first one with spin 1/2, mass, charge and magnetic moment 
Ml, Zi and /xi, respectively, the second one with spin 1, mass, charge, magnetic moment 
and quadrupole moment M2, Z2, 112 and Q2, respectively. The magnetic and the quadrupole 
moments are given in n.m. and fm^, respectively. The non-relativistic reduction of the 
covariant current for the point-like spin-1/2 particle gives for the charge p and current j 
operators in r-space [32]: 

pi(q) = Zie^'i--^ - i ^^l^q • (Si x pOe^'^-^ , 

j^^^^ ^ 2Si^P^'''''"^ " ' ' ^^^^ 

where q is the three-momentum transferred to the particle, pi and Si are the momentum 

and spin operators, respectively, and {•••,•••} denotes the anticommutator. We have here 

neglected the Darwin-Foldy relativistic correction. 

The covariant current operator for a spin-1 particle is written as [33] 

+G2{Q')[e^{e'*-q)-e'*^{e-q)]} , (A2) 
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where E, E' are the initial and final energies, = e{X,pY and e'^ = e(A',p')^ are the four- 
vector spin-1 initial and final polarizations, = p'^ — p^, = p'^ + p'^ and = —q ■ q- 
The three form factors Gi{Q^), G2{Q^) and G3{Q'^) are related to the charge, magnetic and 
quadrupole form factors as [32] 

Gc{Q')^G,(Q') + lrjGQ{Q'), 

Gq{Q') = G,{Q') - G^iQ') + (1 + v)G3{Q') , 

G'm(Q') = G2{Q^) . (A3) 

Here r) = gV(4M|), ^^(O) = Z2, Gm{0) = (M2/M) /i2 and Gq{0) = MfQa, M being the 
nucleon mass. 

To perform the non-relativistic reduction of Eq. (A2), the following relations are used: 
.(A.)-l^,e(A).-HM_,, (A4) 

with e(±l) = ±i, 0), e(0) = (0, 0, 1), and 

1 i 

e(A')ae(A)/3 = 5a/3 - ^{Sa, Sp} + -e„/3^S^ , (A5) 

S being the spin operator. 

The final r-space expressions for the charge and current operators of the spin 1 particle 
are: 



Notations are similar to the ones used in Eq. (Al). It is important to note that besides for 
the quadrupole moment term, Eq. (A6) and Eq. (Al) are formally identical. 

To calculate the MM interaction between the two spin-1/2 and spin-1 particles, we 
consider the standard one-photon exchange Feynman diagram, from which we can write: 
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VMM{r) = j dcie"^"'vMM{ci) , 



vmm{cl) = y^[pi(q)p2(q) - ji(q) • j2(q)] ■ (A7) 

With a straightforward algebra, using Eqs. (Al) and (A6) and keeping terms up to 0(1/M^), 
the formulas for VMM{Nd) of Eqs. (6) and (9) are obtained. 

In an equivalent derivation, VMM{Nd) is written as sum of the NN MM interactions 
between each nucleon of the deuteron and the third particle, at large separation distances. 
It is however important to note that the center of mass (cm.) of each two-body NN 
subsystem is not at rest, and therefore Eqs. (2)-(4), which are derived in the cm. reference 
frame, need to be generalized. In fact, the MM interaction between two spin-1/2 point-like 
particles in a generic reference frame in which the cm. of the system has momentum P, is 
given by [2,34]: 



VMMir) = --Aj^^Su + - - ^)(L . S + L . A) 



+^{-^ ^)(L-S-L-A) 

'-[(rxP)-A+(rxP)-S-^— i]} . (A8) 



4M1M2 M1 + M2 

Here Mj, Zi, /li {i — 1, 2) and M12 are the masses, charges, magnetic moments and reduced 
mass of the two particles, r is their relative position, 5'i2 = 3(Si • r)(S2 • f) — Si • S2 is the 
tensor operator. Si and S2 being the spin operators, S and A are defined as S = Si -|- S2 
and A = Si — S2, L is the orbital angular momentum. The last term of Eq. (A8) is the well 
known Thomas precession (TP) term (see Ref. [35] and references therein). Clearly, Eq. (A8) 
becomes Eqs. (2)-(4), when we consider two nucleons in their cm. reference frame. If the 
TP contribution, which is present only in vmm{pp) {Zi 7^ and Z2 7^ 0), was neglected, the 
p — d MM interaction would have become 
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with same notation as in Eq. (7). 

Finally, note that Eq. (A8) gives the MM interactions also for four-body systems like 
p — ^He and n — ^H. 
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FIG. 1. The n~d Ay calculated using AV18 (solid line) and AV18+MM (dotted-dashed line). 
For the dashed line see text. Experimental points are from Ref. [20]. 
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FIG. 3. The p-d Ay and iTu calculated using AV18 (solid line) and AV18+MM (dotted-dashed 
line). For the dashed line see text. Experimental points are from Ref. [23] (1 MeV) and Ref. [25] 
(3 MeV). 
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FIG. 4. The p-d Ay and iTu calculated using AV18 (solid line) and AV18+MM (dotted-dashed 

line). Experimental points at Ep = 5, 10 MeV and = 10 MeV are from Ref. [27], Ea = 20 MeV 

are from Ref. [28]. 
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FIG. 5. Relative difference between the theoretical and experimental vector analyzing power 
at the peak as a function of energy. 
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FIG. 6. Relative difference between the theoretical and experimental vector analyzing power 
zTii at the peak as a function of energy. 
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FIG. 7. The differential cross section and Ay at Eiab = 65 MeV. For explanation of the curves 
see text. Experimental data are from Ref. [29]. 
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FIG. 8. The deuteron analyzing power iTu and the tensor analyzing powers T2o,T2i,T22 at 
Eiab = 65 MeV. For explanation of the curves see text. Experimental data are from Ref. [30]. 
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